On the Singularity of the Free Energy at First 
Order Phase Transition 



S. Friedli 
Institut de physique theorique 
EPF-L CH-1015 Lausanne Switzerland 

Ch.-Ed. Pfister 
Institut de mathematiques 
EPF-L CH-1015 Lausanne Switzerland 



Abstract: At first order phase transition the free energy does not have an ana- 
lytic continuation in the thermodynamical variable, which is conjugate to an order 
parameter for the transition. This result is proved at low temperature for lattice 
models with finite range interaction and two periodic ground-states, under the only 
condition that they verify Peierls condition. 

1. Introduction 

We study a lattice model with finite state space on Z"^, d >2. Let TYq be a Hamil- 
tonian with finite-range periodic interaction, having two periodic ground-states ipi 
and iIj2, and so that Peierls condition is verified. Let Tii be a Hamiltonian with 
periodic and finite range interaction, so that the perturbed Hamiltonian 

splits the degeneracy of the ground-states of Hq- if /i < 0, then has a unique 
ground-state ip2, and if /i > 0, then has a unique ground-state ipi- The free 
energy of the model with Hamiltonian 7Y^, at inverse temperature (3, is denoted by 
f{fi,P). Our main result is 

Theorem 1.1. Under the above setting, there exist an open interval Uq 9 0, /5* G M"*" 
and, for all (3 > (3* , fi*{(3) G Uq with the following properties. 

(1) There is a first-order phase transition at fi*{/3). 

(2) The free energy f{fi,P) is analytic in n in {n E Uq : fi < it has a 
C°° continuation in {fi E Uq : jj, < fi*{P)}. 

(3) The free energy f{fi,j3) is analytic in fi in {fi E Uq : ji > ji*{[3)}; it has a 
C°° continuation in {fj, E Uq : fi > fj,*{P)}. 

(4) There is no analytic continuation of f from fi < to fi > fi*{l3) across 
fj^*{P), or vice-versa. 



This theorem answers a fundamental theoretical question: does the free energy, 
which is analytic in the region of a single phase, have an analytic continuation 
beyond a first-order phase transition point? The answer is yes for the theory of a 
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simple fluid of van der Waals or for mean-field theories. The analytic continuation 
of the free energy beyond the transition point was interpreted as the free energy of a 
metastable phase. The answer is no for models with finite range interaction, under 



very general conditions, as Theorem p..l| shows. This contrasted behavior has its 



origin in the fact that for models with finite range interaction there is spatial phase 
separation at first order phase transition, contrary to what happens in a mean-field 
model. Theorem |TTT] and its proof confirm the prediction of the droplet model [|I| . 

Theorem |1.1| generalizes the works of Isakov for the Ising model and ||^, where 
a similar theorem is proven under additional assumptions, which are not easy to 



verify in a concrete model. Our version of Theorem which relies uniquely on 
Peierls condition, is therefore a genuine improvement of The first result of this 
kind was proven by Kunz and Souillard |0] ; it concerns the non-analytic behavior of 
the generating function of the cluster size distribution in percolation, which plays 
the role of a free energy in that model. The first statement of Theorem |1.1| is a 
particular case of the theory of Pirogov and Sinai (see 0). We give a proof of this 
result, as far as it concerns the free energy, since we need detailed informations 
about the phase diagram in the complex plane of the parameter /i. 

The obstruction to an analytic continuation of the free energy in the variable fi 
is due to the stability of the droplets of both phases in a neighborhood of fi*. Our 
proof follows for the essential that of Isakov in 0. We give a detailed proof of 
Theorem and do not assume any familiarity with p| or [Q. On the other hand 



we assume that the reader is familiar with the cluster expansion technique. 

The results presented here are true for a much larger class of systems, but for 
the sake of simplicity we restrict our discussion in that paper to the above setting, 
which is already quite general. For example. Theorem |L1] is true for Potts model 
with high number q of components at the first order phase transition point (3c, where 
the q ordered phases coexist with the disordered phase. Here ^ = (3, the inverse 
temperature, and the statement is that the free energy, which is analytic for [3 > (3c, 
or for (3 < (3c, does not have an analytic continuation across Pc- Theorem |1.1| is 
also true when the model has more than two ground-states. For example, for the 
Blume-Capel model, whose Hamiltonian is 

-^(si- Sjf - h^Si- X'^s'^i with Si e {-1,0,1}, 

i,j i i 

the free energy is an analytic function of h and A in the single phase regions. At 
low temperature, at the triple point occurring at /i = and A = \*{(3) there is no 
analytic continuation of the free energy in A, along the path /i = 0, or in the variable 
along the path A = A*. The case of coexistence of more than two phases will be 
treated in a separate paper. 

In the rest of the section we fix the main notations following chapter two of Sinai's 
book IP], so that the reader may easily find more information if necessary. We also 
state Lemma |1.1| which contains all estimates on partition functions or free energies. 
We omit the proof, which relies on the cluster expansion method. 

The model is defined on the lattice l/, d > 2. The spin variables (f{x), x G Z*^, 
take values in a finite state space. If if, ip are two spin configurations, then ip = ip 
(a.s.) means that v?(x) ^ i(>{x) holds only on a finite subset of Z'^. The restriction 
of to a subset A C Z'^ is denoted by (p{A). The cardinality of a subset S is 
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denoted by l^l- li x,y G Z'^, then \x — y\ := maxf^^ |xj — y^l; ii W d 1/ and 
X G Z'^, then 14^) := min^g^ |x — y\ and if W,W' are subsets of Z*^, then 
d{W, W) = min^6H/ W). We define for C Z'^ 

dW ■.= {xeW : d{x, Z'^\W) = 1} . 

A subset W C Z'^ is connected if any two points x,y E W are connected by a path 
{xq, Xi, . . . , Xn} C W, with Xq = X, Xn = y and |xj — Xi+i\ = 1, i = 0, 1, . . . , n — 1. 
A component is a maximally connected subset. 

Let be a Hamiltonian with finite-range and periodic bounded interaction. By 
introducing an equivalent model on a sublattice, with a larger state space, we can 
assume that the model is translation invariant with interaction between neighboring 
spins ip{x) and ip{y), \x — y\ = 1, only. Therefore, without restricting the general- 
ity, we assume that this is the case and that the interaction is Z''-invariant. The 
Hamiltonian is written 

= no + fiHi , /i G M . 

Tio has two Z'^-invariant ground-states ipi and ip2, and the perturbation Tii splits 
the degeneracy of the ground-states of Hq. We assume that the energy (per unit 
spin) of the ground-states of TCq is 0. U!^{(p) = Wo.x + /^^i.a: is the interaction energy 
of the spin located at x for the configuration if, so that by definition 

= ^xi^) (formal sum) . 

Ui,x is an order parameter for the phase transition. If and ijj are two configurations 
and if = ip (a.s.), then 

This last sum is finite since only finitely many terms are non-zero. The main condi- 
tion, which we impose on Tio, is Peierls condition for the ground-states ipi and il)2- 
Let X G Z^ and 

Wi{x) := {yeZ'^: \y-x\<l}. 
The boundary dip of the configuration ip is the subset of Z"' defined by 

dip:= [j {W^{x)■. y^{W^{x)) ^ tPmiWiix)) , m = l,2}. 

Peierls condition means that there exists a positive constant p such that for m = 1, 2 
'Ho{ip\ipm) > pldipl V (p such that (p = tpm (a.s.) . 

We shall not write usually the /^-dependence of some quantity; we write for example 
H or Ux instead of or W^. 

Definition 1.1. Let M denote a finite connected subset ofl/', and let ip be a config- 
uration. Then a couple T = (M, ip{M)) is called a contour of (p if M is a component 
of the boundary dip of ip. A couple T = (M, ip{M)) of this type is called a contour if 
there exists at least one configuration ip such that T is a contour of ip. 

If r = (M, ip{M)) is a contour, then M is the support of F, which we also denote 
by suppF. Suppose that F = (M, y9(M)) is a contour and consider the components 
Aa of 'lf-\M. Then for each component Aq, there exists a unique ground-state 
'^q{a)i such that for each x G dA^ one has ip{W\{xy) = ipq{a)(Wi{x)). The index 
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q{a) is the label of the component A^,. For any contour F there exists a unique 
infinite component of Z'^ysupp F, Ext F, called the exterior of F; all other components 
are called internal components of F. The ground-state corresponding to the label 
of Ext F is the boundary condition of F; the superscript q in F'' indicates that 
F is a contour with boundary condition ipq. Intm F is the union of all internal 
components of F with label m; Int F := lJm=i 2 I^^m F is the interior of F. We use 
the abbreviations |F| := |supp F| and Kn(F) := |Int„F|. We defineQ 

(1.1) V{T''):=Vm{T^) m^q. 
For X e Z'^, let 

c(x) := {y : max \xi - y.i\ < 1/2} 

i=l 

be the unit cube of center x in M"^. If A C Z*^, then |A| is equal to the (i- volume of 

(1.2) |Jc(x)cM^ 

The [d — l)-volume of the boundary of the set (|1.2| ) is denoted by d\h\. We have 



(1.3) 2c/|A|— <9|A|. 

The equality in ( |1.3| ) is true for cubes only. When A = Int^F'', m 7^ g, V{V^) = |A| 
and dV{r'^) = <9|A|; there exists a positive constant Co such that 

(1.4) dV{r'^) < Col^l q=l,2. 

For each contour F = {M,ip{M)) there corresponds a unique configuration ifr with 
the properties: (fr = "^q on ExtF, where q is the label of Ext F, ifr{M) = if{M), 
V^r = on Intm F, m = 1, 2. F is the only contour of ipr- Let A C Z'^; the notation 
F C A means that suppF C A, Int F C A and (i(supp F, A'^) > 1. A contour F of a 
configuration cp is an external contour of (/? if and only if F C Ext F' for any contour 
F' of ^. 

Definition 1.2. Let f2(F'^) be the set of configurations ip = ipq (a.s.) such that 
is the only external contour of (f. Then 

e{T^):= Yl exp[-/?H(^|^,)]. 

Let K d J/- he a finite subset; let fig (A) he the set of configurations (f = ipg (a.s.) 
such that F C A whenever T is a contour of ip. Then 

e,(A):= exp[-/57i:(^|^,)]. 
Two fundamental identities relate the partition functions 6(F^) and Qq{A). 

n 

(1-5) 0,(A) = En0(rn, 

i=l 

where the sum is over the set of all families {F^, . . . , F^} of external contours in A, 
and 

2 

(1.6) e(F'^) = exp [ - pni^prAM H ©^(Int^F^) . 



m=l 



^Here our convention differs from 



We define (limit in the sense of van Hove) 

gq := lim --r^ log Qq{k) . 

ATZ'* I3\A\ 

The energy (per unit volume) of ifjm for the Hamiltonian Tii is 

By definition of Tii, /i(^/'2) ~ ^(^i) 7^ 0, and we assume that 

A := h{^2)-h{^i) > 0. 
The free energy in the thermodynamical limit is 

(1.7) f = lirn --l-\oge,{A) + lim -^Y^U.ii^,) = 

' ^' ' III 

It is independent of the boundary condition ijjq. 

Definition 1.3. Let be a contour with boundary condition ipq. The weight uj(V'^) 
of V is 

.(n:^exp[-,«(.„K)] n 

The (bare ) surface energy of a contour is 

For a contour we set 

xGsupp ri 

Since the interaction is bounded, there exists a constant Ci so that 

(1.8) |a(y.rOI <Ci|r''|. 
Using these notations we have 

xGsupp ri xeint r'' 

= Hoi^r^l^g) + /ia(<^rO + KHiPm) - /i(^,))K(n) 

(1.9) =||n||+/ia(v.rO+MM^J-MV'g))V^(r'') (m^g). 

The surface energy Hr*^!! is always strictly positive since Peierls condition holds, and 
there exists a constant C2, independent of g = 1, 2, such that 

(1.10) p|r«| < lir^ll < Csl^l . 

Definition 1.4. The weight ujij"^) is r-stable for V if 

\uiT'')\ < exp(-r|r^|) . 

For finite subset A C Z*^, using ( |1.5| ) and ( |1.6| ), one obtains easily the following 
identity for the partition function Og(A), 



1), 



(1.11) e,(A) = i + 5^n^(n 

i=l 

where the sum is over all families of compatible contours {F^, . . . , F^} with boundary 
condition ipq, that is, F^ C A and (i(supp F^, supp Fp > 1 for alH 7^ j, 2, j = 1, . . . , n. 



n > 1. If the weights of all contours with boundary condition ipg are r-stable and 
if T is large enough, then one can express the logarithm of 6g(A) as an absolutely 
convergent sum, 



^ m 

[1.12) ioge,(A) = 5^-5^--- ^^(r?,...,r^)n^(rr 



m>l r«cA r«,cA i=l 



In ( |1.12|) (pJ^{Tl, . . . , r^) is a purely combinatorial factor. This is the basic formula 
which is used for controlling Qg{A). We also introduce restricted partition functions 
and free energies. For each n = 0, 1, . . . , we define new weights ujn(r'^) 



uj{T'i) if ViV^) < n, 
otherwise. 



:i.l3) ^«:=_lim^loge^(A) and f^:=g^ + zh{%). 



For q = 1,2, we define 6^ by equation replacing uj(T'^) by ci;„(r^), and we 

set (provided that 6^ (A) ^ for all A) 

1 

atP i3\A\ 

fg is the restricted free energy of order n and boundary condition ijjq. Let 
(1.14) l{n) := Cq^ \2dn^] n>l. 

Notice that e^(A) = Qq{A) if |A| < n, and that 1/(r«) > n implies that |r«| > l{n) 
since (O) and (ITl hold. 



Lemma 1.1. Suppose that the weights (^(F'') are r-stable for all V^. Then there 
exists Kq < oo and Tq < oo, so that for all t > Tq, ( |1.12|) is absolutely convergent 
and 

P\gg\<Koe--. 

For all subsets A C Z^, 

\logQg{A)+/3gg\A\\<Koe~^d\A\. 
//cj(F«) = for all such that |F«| < m, then 

f3\9,\ < {Koe-^r. 

For n > 1 and m > n, 

Furthermore, if uj{T'^) depends on a parameter t and 

|^^(F'')| <Die--|r^l and | ^c.(F'') | < D^e'-I^'l , 

d^ 

then there exists < oo and < oo, k = 1,2, so that for all t > r^, ~^9q exists 
and 

P\j^g,\ < D^Kie-^ and /5| < max{D2, D^jA'ae^^ 

For all subsets A C 1/ , 

1 1 log e,(A) + [5j^gg I A| I < D,K,e-^ d\A\ 
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and 

I ^ log e,(A) + [3^g, I A| I < max{D2, D^jK^e-^ d\A\ . 



Lemma is proved by the cluster expansion method. It follows from (|1.12| ) and 



arguments similar to those of the proof of Lemma 3.5. in section 3.3 in ||^. 



2. Proof of Theorem 1.1 



The proof of Theorem |1.1| is given in the next five subsections. In subsection 



2A\ we construct the phase diagram and in subsection p.2| we study the analytic 
continuation of the weights of contours in a neighborhood of the point of phase 
coexistence /i*. The results about the analytic continuation are crucial for the rest 
of the analysis. Construction of the phase diagram in the complex plane has been 
done by Isakov We follow partly this reference and Zahradnik [0. In subsection 
|2.3| we derive an expression of the derivatives of the free energy at finite volume. 
We prove a lower bound for a restricted class of terms of this expression. This is an 
improved version of a similar analysis of Isakov [Q] . From these results we obtain a 
lower bound for the derivatives of the free energy /a in a finite box A. We show in 
subsection p.4| that for large (], there exists an increasing diverging sequence {kn} 



d 



SO that the /c^'^-derivative of /a with respect to /i, evaluated at behaves as 
(provided that A is large enough). In the last subsection we end the proof of the 
impossibility of an analytic continuation of the free energy across /i*, by showing 
that the results of subsection p.4| remain true in the thermodynamical limit. 



2.1. Construction of the phase diagram in the complex plane. We construct 
the phase diagram for complex values of the parameter /i, by constructing iteratively 
the phase diagram for the restricted free energies (see ( p..l3D ). We set 2; := n + iu. 
The method consists in finding a sequence of intervals for each 1/ G M, 

Ur.{v- /3) := ifiUu; (3) - hi fi^iu; /3) + bl) , 

with the properties 

(2.1) (/i:(z/; P) - bl filiu- P) + bl) C ifiUiu; /5) - Ci, /i:_i(z/; P) + Ci) 

and lim^fe^ = 0, g = 1,2. By construction of the intervals f/„_i(z/;/5) the restricted 
free energies /^~^ of order n — 1, q = 1,2, are well-defined and analytic on 

U„_i ■.= {z eC: Rez e f/„_i(Imz; /3)} . 

The point /i* (z/; /3), n > 1, is solution of the equation 

Re{fr'mu; P) + iv) - /r-i(/i:(z/; P) + iv)) = . 

yu* (0; P) is the point of phase coexistence for the restricted free energies of order n—1, 
and the point of phase coexistence of the model is given by /i*(0; P) = lim„ //* (0; P). 
This iterative procedure also gives the necessary results needed in subsection |2.1| 
about the analytic continuation of the weights ujIV^) around the point of phase 
coexistence /i*. Since we need sharp results about the analytic continuation of the 
weights ci;(r''), we must choose carefully the two sequences {&^}, q = 1,2. In order 
to ease the exposition we first describe the iterative procedure with a specific choice 
of {6^}, based on the isoperimetric inequality 

(2.2) ViT")'^ <x^^\\T''\\ Vr«,g=l,2. 
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Existence of x ( P-2| ) follows from ( |1.3| ), ( |1.4| ) and ( |1.10| ). Then, in subsection |2.2| , 



we make another choice for {fc^}. This iterative construction is given in details in 



the proof of the Proposition which is the main result of subsection 2.1 . 
Proposition 2.1. Let < e < p and < 5 < 1 so that A - 25 > 0. Set 
Uo ■= {-C^^e, C^h) and Uq := {z G C : Re^ G Uq} 

and 

t{(3) ■.= f3{p-e)-3Co5. 

There exists /3o G IR"*" such that for all (3 > [3q the following holds. 

A) There exists a continuous real-valued function on R, z/ i— > /i*(z/;/5), so that 
fi*{iy;p) + iu G Uq. 

B) If fi + iu G Uo and p < i/ien the weight ci;(r^) is T{l3)-stable for all 
contours with boundary condition ip2, and analytic in z = p + iu if ji < p*{iy; P). 

C) If p + IP E Uo and p > p*^^;^), then the weight uj(T^) is T{j3)-stable for all 
contours with boundary condition ipi, and analytic in z = p + iv if p> p*{u\ (3). 

Remark 2.1. p is the constant of the Peierls condition and A = /i(V^2) ~ ^(V^i) > 0. 
We may choose 5 in such a way that 5 = 6{P) and lim^^oo ^iP) = 0, without 
changing the theorem. Indeed, the only condition which we need to satisfy is ( |2.6| ). 
So, whenever we need it, we consider 6 as function of /5, so that by taking j3 large 
enough, we have 6 as small as we wish. 

Proof. The iterative method depends on a free parameter 6'', < < 1. On the 
interval Uq^u; P) := {—bo, bo) with bo = eC^-^, fg{p + iu) is defined and we set 
Poii^; P) := 0. The two decreasing sequences {fc^}, q = 1,2 and n > 1, are defined in 
(p.8|). The iterative construction is possible whenever the sequences {6^}, q = 1,2, 
verify ( |2.7| ), ( |2.13| ) and ( |2.14| ). We prove iteratively the following statements. 

A. fg{p + iu; P) is defined for all p G Un~i{u; P), and u i— > /i*(z/; P) is a contin- 
uous solution of the equation 

Re(/r^(/i:(z/; p) + lu) - fr\i^n{^; p) + ^^)) = o , 

so that holds. 

B. On U„, co'„(r^) is analytic for any contour V^, g = 1,2, and Un{T'^) is ti(/3)- 
stable (see ( p.3|) ). 

C. On U„, |-?^cj„(n)| < /5C3e-"^(^)l^'l (see (O) and (|D). 

dz 

D. For each n > 1, ii p < p*^{u;P) — b\, then uj{T'^) is r(/3)-stable for any F^ 
with boundary condition ip2. Similarly, for each n > 1, ii p > Pn{u; P) + b\, 
then uj{T^) is r(/5)-stable for any F^ with boundary condition ipi. 

From these results the proposition follows with 

l^*{^]P) = lim Pn{jy;P) ■ 

n—>oo 

The analyticity of the weights co'(F'') is an immediate consequence of their stability 
since 6m,(IntmF^) and 9q(IntmF'') 7^ are analytic. 

Let < 6'' < 1 be given, as well as e and 6 as in the proposition. We introduce 
all constants used in the proof below. 

(2.3) niP; e') := P{p{l - 6') - e) - 26Co , 
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d 



(2.4) 
and 

(2.5) Cs := Ci + 26Co + (A + 25)ix~^C2)^^ . 

We assume that /?o is large enough so that|r2(/3) > max{r*,r*,r*}, (pl8|) holds, 

(2.6) i^e-"i(^)<5 and Cgi^e-'^^^) < 5 , 



where -fC = max{KQ, Ki}, and i^o? -^i are the constants of Lemma |1 . 1| . We assume 
that for q = 1,2, 

2X/{n) 

(2.7) - - 



If we define 
(2.8) 



61 = hi ■-- 



\/n> I. 



n > 1 



(A + 2d)n^ 

then it is immediate to verify { \2.7\ ) when (3 is large enough or 6 small enough. On 
Uq all contours F with empty interior are f3{p — £:)-stable (see ( p..8[ )), and 

^^(r) < /3Ci|r|e-^(^-^)iri < /5Cie-[^(''-^)-^liri < pC,e-^'^^^\^\ . 
dz 



Assume that the construction has been done for all m < n — 1. By Lemma |LT 
if 2; G U„_i, then 



(2.9) 



-^qT\ < C,i^e-"^(^) <6 m<n-l. 



' dz 



A. We prove the existence of /i* (i^; P) G U„„i. /i* (z^; P) is solution of the equation 
Let F'^iz) := f^{z) - f^{z). Then, for fi' + ii^ E U„_i, 



(2.10) 



:^n—l/,/ 



{jj.' + if) 



?n-l/- ,,/ 



(/i' + iu) 



(/i*_i + ijy) 



F""i(/i;_i + zz/) 



= F"-i(^' + zz/) - F""i(/i:_i + lu 

r' A 

If V^(r) = 77,-1, then |r| > /(ra — 1). Therefore, by Lemma [1.1| , 

(2.11) |(^n-l_^n-2)(^*_^^^^)| 

If z' = /i' + zi^ G U„_i, then ( p:TD| ), (|2;|) and ( pi]) imply 
A(/x' - /i:_i) + 25|/i' - + 2/5-i5'("-i) > ReF"-i(2') 

> A(/i' - - 25|/i' - - 2/3-i5'("-i) 



^Tj!, /c = 0, 1, 2, are defined in Lemma |l . l| . Condition T2(/5) > is needed only in L emma 2.2 



We have stated Lemma 2.2 separately in order to simplify the proof of Proposition 2.1 
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( p.TD implies 



^ —W~> 



2^Z(n-l) 



/3(A - 25) ' 

SO that ReF"-i(/i;_i - b^^^ + iu) < and ReF"-i + bl_^ + iz/) > 0. This 
proves the existence of yU* and its uniqueness, since fi ^ ReF"~^(/i + iu) is strictly 
increasing. Moreover, by putting /i' = //* (i^; /3) in ( |2.1CI| ), we get 

2^/(n-l) 

Therefore U„ C U„_i. The implicit function theorem implies that u i— > /z* (i^; /3) is 
continuous (even C°°). 



B. We prove that u;„(r'') is ri-stable for all contours V^, q = 1,2. Let be a contour 
with VlV^) = n. All contours contributing to 0m(lntmr^) and 0g(lntmr'') have 
volumes smaller than n — 1, so that for these contours ci;(r) = ct;„_i(r). If 2; G Un-i 
(use ([T7^ ), ( |1.4D and the definition of [/q), then 

(2.12) 

e„(int„n) 



\L0 



(n)|=exp[-/3ReH(yPr.|^,)] \{ 



e,(int™r'?) 



< exp 
= exp 



ir^ll + {f3e + 2Co6) \T^\ - /3Re(/r' - /r')^(r' 

v(r'^) 



Let /i G Un-i{i^',P)- We prove that 6^ verify the following conditions, which imply 
the Ti-stability. 



(2.13) 



(2.14) 



-Re{fr' - /2 -')^ <0' if /X < + 6^ and ^(r^) 



-Re(/. 



n-1 



/r-^) ^ <9' iffi>fi:-bl and r(ri) 



n, 



n. 



For the present choice of {b^}, the isoperimetric inequality (|2.2|) implies 



X 



and therefore 

|Re(/r^-/r^) 



,1/(1") 


Re f 


lir^ii 


J n 



<l/i-/i;i(A + 25)^< 



ir^i 



Conditions ( p.l3| ) and (|2.14]) ensure that on U„ 

|cj(n)| <exp [-/3(p(l-^')-^-2/?"^Co(5)|n 



C. We prove that on I[J„ 

|-^...(r)|</3C3e--(^)l^l. 
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Let 1/(19) = n; from (|T;| 
d 



dz 



^(ioge„(int„n) -ioge,(int^n))) . 



By Lemma IJ. and the isoperimetric inequality ( |1.3|) we get 



(2.15) 



< f3\ujnin\{\mc^ + 26Co) + VimA + 26)) 

< /?C3e-^^(^)ir'l . 



D. We prove that uj{r'^){z) is r(/5)-stable for any contour with boundary condition 
4'2, if < f^ni^'^ P) ~ ^n- Using the induction hypothesis it is sufficient to prove this 
statement for z = /i + zz/ G U„_i and /i < fJ^n^i^] (3) — h\. \i z = ^ + iv E U„,-_i, 
then all contours with volume V{T) < n — 1 are ri(/3)-stable, and for < /i*, 
/i I— > Re(/"~^ — f2~^){l^ + ^'^) is strictly decreasing. If /i < /i^(z^; /3) — then (see 
O and (^)) 



/rO(/i + zi.) = -/3 



d^ 



Re(/r^-/r')(/i+2j/)c?/x 



(2.16) 



> -P / ;^Re(/r^ - /r^)(/i + lu) d^i 

> phl{A - 25) > 25'(") . 



First suppose that V{T^) < n. From (^?I6D and (|]T2]) it follows that u;{T^) is 
/3(p — £ — 2/3~^Co5)-stable, in particular r(/3)-stable. Moreover, if |A| < n, then 



(2.17) 



exp[-(3z{h{^|J,)-h{^|J2))\A\] 



0i(A) 



e2(A) 



< e' 



,3<5a|A| 



Indeed, all contours inside A are ri(/3)-stable. By Lemma |TTT 

-/3^(fe(V>i)-fe(t/.2))|A| Ql(A-) 



e2(A) 



g-/3Rc(/i"-i(2)-/2"-i(2))|A|g25S|A| 



< e 



2(59|A| 



To prove point D, we prove by induction on |A| that ( [^.17] ) holds for any A. Indeed, 
if (|2.17| ) is true and if we set A := IntiF^, then it follows easily from the definition 
of cu(F2) and from (g^) that ^^(F^) is r(/3)-stable. 



The argument to prove ( p.l7| ) is due to Zahradnik |^. The statement is true for 
|A| < n. Suppose that it is true for |A| < m, m > n, and let |A| = m + 1. The 
induction hypothesis implies that uj{r'^){z) is r(/3)-stable if V^T"^) < m. Therefore 



-/32(h(V>i)-MV'2))|A 



I 01(A) 
e2(A) 



^ \^-f3(zh(i,,)-zh{^2)-9T)\MQ^(^/^^Ud\A\ _ 
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From ([LSP 

0i(^) = Eri0(r])> 

where the sum is over all families {Fj^, . . . , Fj} of compatible external contours in A. 
We say that an external contour Fj is large if ^(Fj) > n. Suppose that the contours 
F]^, . . . Fp are large and all other contours F^^^, . . . Fj not large. We set 

p 

Ext?(A) := (p|ExtF]) n A. 



Summing over all contours which are not large, and using (|1.6| ), we get 



0i(A) = ^eri(Ext?(A)) J]exp[-/37^(<^r]l^i)]0i(IntiF])e2(Int2F]) 
ei(IntiFi) 

.^ij^e.(M,r;)e,(i„t.r;); 

the sums are over all families {Fj, . . . F^} of compatible external large contours in A. 
All contours which are not large are ri(/3)-stable, and we use the cluster expansion 
to control erHExt?(A)), GallntiF]) and BslIntaF]). Notice that (9|Ext?(A)| < 
9|A| + Xlj=i C*o|r]|- By Lemma |TTl| and the induction hypothesis, 

e2(A)|-' 

. J^g-(/3p-/3e+6C7o<5)|rj|g-/3Re(/f-i-/2"-i-32"+92""')|rjl _ 

i=i 

We define 

f{(3) ■.= P{p-e)-6CoS. 
From ( |2.16| ) and Lemma |L1| we have 

- /r' - 9T + gr') > ^'^'^'^ ■ 

Hence, 



g-/32(h(V'i)-M'/'2))|A 



0i(A) 



e2(A) 

We define 



^ g259|A| ^g-5Hn)|ExtP{A)| -Qg-(5'(")+f(/3))|ri| 
g-(f(/3)-Co5)|r| if |r| > 

otherwise. 



Let 0(A) be defined by ( |1.11| ), replacing ti;(F'^) by cj(F), and let 

ATZ'* p|A| 

We assume that jSg is large enough so that for all P > Po, 
(2.18) Ke-^^f^^<6, 
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where K is the constant of Lemma |L^. Since (3\g\ < putting into evidence a 
factor e^^'"^' , we get 

p 



(2.19) 



-pz{h{i,^)-h{^2)m 



0i(A) 



e2(A) 



p 



r(/3)|riL -/3g|Intri| 



In the last hne of ( p.l9| ) we interpret e"^^'^"^**"^' as a partition function (up to a 



boundary term), since by Lemma |LT 



e-/3§|intri| <0(intri)e^»'^|r'l 
We sum over external contours in ( p.igp and get 



-(3z{h(:^i)-h{i,2))\h-\ 



0i(A) 



e2(A) 



< e 



2<5a|A|+/3glA| 



e(A) < e' 



M9|A1 



□ 



It is not difficult to prove more regularity for the curve v t— > /i*(z/;/5). We need 
below only the following result. 

Lemma 2.1. Let < 5 < 1. If 13 is sufficiently large, then for all n > 1 
— /x;(0;/3) = 0, and 



|^M-;/5)|< 



25 \ 2 26 



1 . 



A-26\\A-26J A -26 

Proof. Let 6 be as in the proof of Proposition Because the free energies 
and /^~^ are real on the real axis, it follows that u k-s- /i* (z/; /5) is even, and therefore 

— /X* (0; P) = 0. By definition //* (z/; /5) is solution of 
du 

Re(/r^(/x:(i^; /3) + iv) - fr\f^n{^: /3) + ^i^)) = , 
which implies that 



dv dfi 



and 



- ^ + |^Re«- - gV) . 



d^dv dv du 

From the proof of Proposition |2.1| we have on Um, 

, d 



Let T3(/9) := Ti(/9) — 2^3Y- A similar proof shows that for P sufficiently large, there 
exists C4 so that for any m 

|^^^(r)|</5^C4e--(^)'^l. 
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Assume that /5 is large enough so that 

/3max{C4,C|}i^2e-"^(^^l^l <6. 
Then Lemma |ri| gives for := Re(^['~^ - g'^'^) 

j2 j2 j2 

' < 25 , I— < 25 , l-^G"^^! < 25. 



Hence 

25 // 25 Y 25 
' rfi/ ' - A - 25 ' I I - A - 25 UA - 25/ A - 25 

□ 

Proposition 2.2. Under the conditions of Proposition t/iere exist /5o G IR"^ 
p G N so i/iai i/ie following holds for all (3 > (3o- Let 

r'iP) :=r(/5)-max{^,p}. 

A) If fi + iu ^ Vq and /i < fi*{iy; 13), then 

|^a;(n(.)|</5C3e-'(^)'^\ 

B) If ^ + iu E Uo anc? > P), then 

|-^a;r)(.)|</5C3e--'(^)l^r 

Proof. We consider the iterative construction of the proof of Proposition with 
the same choice of the sequences {6^}. Suppose that z = ^ + iu E U„_i\U„ and 
/i < p). Suppose that ViV"^) < n. We get (see (|2J[5|) ) 

' dz 

Since by Proposition |2.1| u;(r^) is r(/3)-stable, we get for all such that V^(r^) < n, 

|^^(r2)| < /3C3|r2|^e--('^)ir'i < /?C3e--'(^)ir'i . 

Suppose that ViV"^) > n+1. We estimate the derivative at z of ci;(r^) using Cauchy's 
formula with a circle of center z contained in {/i + zz/ : /i < /?)}. We estimate 
from below |Rez — when 2; G U„_i\U„, uniformly in v. 

\Rez - > |Rez - /i^l - l/i^ - > fe^ - l/i^ - . 

We estimate — by first estimating — /i* |. Let m > n; then, since G U„, 

= Re{fr\f^*j - fr\f^*j) - Mfr'ii^:) - fr'if^:)) 
= Re{fr\f^*j - fr\f^*j) - Re{fr\f^*j - fr\f^*j) 

From (|2ll| ) we get 
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SO that 

(2.20) |/i*(z/;/?)-/i:(i.;/3)|< 



/5(A - 26) ■ 

If V{T^) > n + 1, then \T^\ > l{n + 1). Choose p G N so that for all n > 1 

< — ^ <—^^--^-lTn^^<bl-\f^*-f^:\<\Rez-fi* 



|r2|P - V2rfn^ ^ - (A + 26)n-^ /^(A - 25) 
We use Cauchy's formula with a circle of center z and radius |r^|~P and get 

' dz ' 



□ 



2.2. Analytic continuation of the weights of contours at /i*. In this subsec- 
tion we consider how the weight a;(r^) for a contour with boundary condition ip2 
behaves as function of z = + w in the vicinity of z* := /3) + iu. We improve 
the domains of analyticity of the weights of contours, by making a new choice of the 
sequences {fc^}, q = 1,2. The main result of this subsection is Proposition ^]B|. At 
z* the (complex) free energies fq, q = 1,2, are well-defined and can be computed by 
the cluster expansion method. Moreover, 

Re/2(^*) =Re/i(^*). 

Therefore 

Reg,{z*) + /i*(i/; /3)/i(V'i) = Reg,{z*) + /i*(z/; /?)/i(^2) • 
With 6 as in the proof of Proposition we get 



l/i*(^;/5)l<||, 

and 

\u;{T^){z*)\ < exp [ - P\\T^ + + 6Co\T'^\] , V P^ 

We set 

/i*:=/i*(0;/3), 

and adopt the following convention: if a quantity, say H or /g, is evaluated at the 
transition point /i*, we simply write Ti* or /*. 

The analyticity properties of co'(P^) near /i* are controlled by isoperimetric in- 
equalities 

(2.21) V{T^)^ <X2{n)~^\\T^\\ yT\V{T^)>n. 

The difference with ( |2.2| ) is that only contours with boundary condition and 
^(r^) > n are considered for a given n. By definition the isoperimetric constants 
X2{n) verify 

X2{n) 



^ := inf : ^^^^j^^" < C , V P^ such that V{T^) > n| . 

X2('^) is a bounded increasing sequence; we set ^2(00) := lim„X2('^)5 and define 

X2im) 



i?2(w) := inf 



m:m<n 



md 
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There are similar definitions for Xii''^) ^"^^ Ri{n). Tlie corresponding isoperimetric 
inequalities control the analyticity properties of u;(r^) around fi*. 

Lemma 2.2. For any x'q < Xg(oo)j there exists N{x'g) such that for all n > N{x'q), 



For q = 1,2, nRq{n) is increasing in n, provided that a > ^■ 

Proof. Let q = 2 and suppose that 

R2\n) = 1 — for m < n. 

rrid 

Then R2{m') = R2{n) for all m < m' < n. Let n' be the largest n > m such that 
We have n' < oo, otherwise 



1 

nid 



< R2{m) = i?2(n) < Wn>m 



which is impossible. Therefore, either 

= ^ or fl,(„' + i) = 2La(!!:±il, 

n' d {n' + l)d 

and for all A; > n' + 1, since X2{'m) is increasing, 

(2.22) i?,(fc)=inf^= inf ^> inf ^^(^') - ^^(-') 



m<k mil n'<m<k j^^S n'<m<k jji^ 



Inequality (|2.22|) is true for infinitely many n'; since there exists m such that X2 ^ 



X2(^), the first statement is proved. 

On an interval of constancy of R2{n), n 1— > n°'R2{n) is increasing. On the other 
hand, if on [mi,m2] 

X2{n) 



R.in) 



1 1 
rid 

1 



then n 1— n"'R2{n) is increasing on [mi, 7712] since n t— > X2(^) and n n"" d are 
increasing. 

□ 

The next proposition gives the domains of analyticity and the stability properties 
of the weights uj{T) needed for estimating the derivatives of the free energy. 

Proposition 2.3. Let < 9 < 1 and < £ < 1 so that p{l - 9) - e > 0. 

There exist 0<6<1,0<9'<1 and f3o G IR+, such that for all (3 > (3o co'(r^) is 
analytic and ti{/3; 9')-stable in a complex neighborhood of 

{zeC: Rez< i2*{lmz; (3) + 9A^^R2{V{T^))} CMJo . 

Moreover 

|-^o;(r2)| </3C3e--^(^^^')|rr 
' dz ' 
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Similar properties hold for u;(r^) in a complex neighborhood of 

{zeC: n*{lmz; (3) - eA-^Ri{V{T^)) < Rez} f] Uq . 



Ti{(3;6') and T2{(3;6') are defined at (U^ ) and (|2] 



Proof. If in the iterative method of the proof of Proposition |2.1| we find < 9' < 1 
and bl, bl, so that ( CT) and 

(2.23) {fi*iiy;f3)-eA~'Ri{n),fi*iu;l3) + eA-'R2in)) Cf/„(i/;/5) 

hold, then Proposition ^]3| is true. Formula (|2.23|) is satisfied if (see ( |2.2CI|) ) 



bl > eA~^R,{n) 



and this is the case if 



/3(A - 25) 
C 



bl := QA-'R,{n) + -S 



with C a suitable constant, which is chosen so that ( |2.7| ) is also satisfied. If P is 
large enough and S small enough, then there exists ^' < 1 so that ( |2.13|) and ( p.l4|) 
hold. Indeed, let V^(r^) = n, z = fi + iu and /i < /U.*(i/; P) + b\] then 



ir^ll X2{n) 

- A P X2{n) - 



□ 



2.3. Derivatives of the free energy at finite volume. Although non-analytic 
behavior of the free energy occurs only in the thermodynamical limit, most of the 
analysis is done at finite volume. We write 

t=t' 

for the k^^ order derivative at t' of the function g. The method of Isakov allows 
to get estimates of the derivatives of the free energy at /x*, which are uniform in 
the volume. We consider the case of the boundary condition '02- The other case is 
similar. We tacitly assume that (3 is large enough so that Lemma |TTT| and all results 
of subsections |2.1| and ^.2| are valid. The main tool for estimating the derivatives of 
the free energy is Cauchy's formula. However, we need to establish several results 
before we can obtain the desired estimates on the derivatives of the free energy. The 
preparatory work is done in this subsection, which is divided into three subsections. 



In p.3.1| we give an expression of the derivatives of the free energy in terms of the 
derivatives of a free energy of a contour M(r^) = — log(l + 0A(r^)) ^ — 0A(r^) (see 
( p.25| )). The main work is to estimate 

The boundary of the disc Dr is decomposed naturally into two parts, dD^ and dD!^, 

and the integral into two integrals Il^{T^) and I^JT^) (see and (^^)). In 

d r-r2\ 



2.3.2| we prove the upper bound ( |2.29D for /|^(r^), and in |2.3.3| we evaluate /^^(r 
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by the stationary phase method, see ( |2.33| ) and ( p.34| ). This is a key point in the 
proof of Theorem since we obtain lower and upper bounds for /^^(F^). 

2.3.1. An expression for the derivatives of the free energy. Let A = A(L) be the 
cubic box 

A(L) := {zeZ'^ : \x\ < L} . 

We introduce a hnear order, denoted by <, among all contours F*^ C A with bound- 
ary condition ipq. We assume that the linear order is such that V{r''^) < l^(F'^) if 
Y'<i < T''. There exists a natural enumeration of the contours by the positive inte- 
gers. The predecessor of F'' in that enumeration (if F"^ is not the smallest contour) 
is denoted by i(T'^). We introduce the restricted partition function Gr9(A), which 
is computed with the contours of 

Ca(F«) := {F"? C A : F'" < F"} , 

that is 



1), 



(2.24) er,(A):=l + 5^n^(n 

i=l 

where the sum is over all families of compatible contours {Ff , . . . , F^} which belong 
to Ca(F^). The partition function 6q(A) is written as a finite product 

By convention 6j(r'3)(A) := 1 when F'' is the smallest contour. We set 

Wi(r9)(Aj 

u\(r'^) is the free energy cost for introducing the new contour F*^ in the restricted 
model, where all contours verify F"^ < F'^. We have the identity 

er,(A) = e,(r,)(A) + uiT") e,(r,)(A(F'')) 

0.(r.)(A(F'?))- 



= e,(r.)(A) 1+.;(F'^) 

V fc)i(r9)(A) 

In this last expression 0j(rg)(A(F'^)) denotes the restricted partition function 



e,(r,)(A(F«)):=l + ^n^(r: 



1=1 



where the sum is over all families of compatible contours {Ff , . . . , F^} which belong 
to C\{i{T'^)), and such that {F^, Tf, . . . , F^} is a compatible family. We also set 

With these notations 

(—1)"' 

(2.25) u^{T^) = - log (1 + Mn) = E -^M^"" , 

n>l 

and for k > 2 

TIC A 
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We consider the case of the boundary condition '?/^2- [0A(r^)"]^*«'' is computed using 
Cauchy's formula, 

where dD^. is the boundary of a disc of radius r and center /i* inside the analyticity 
region of Proposition |2.3|, 



UoH G C : Rez < /i*(Im(z); /5) + M-^i^slV'lr^))} 
The function real on the real axis, so that 



0A(r2)"(^) _ 0A(r2)"(^ 

and consequently 



Remark 2.2. From Lemma |2]l|, there exists C" independent of v and n, so that 

/i:(^;/5)>/i:(0;/3)-CV. 

This implies that the region {Rez < fi* — C"(Im2;)^ + 6'A~^i?2(V^(r^))} is always in 
the analyticity region of w(r^), which is given in Proposition |2.3| . Therefore, if 



C'< 



2{9A-^R2{V{T^))) 



2 ' 



then the disc Dr of center /i* and radius r = 6'A ^i?2(^(r^)) is inside the analyticity 
region of uj{r'^). This happens as soon as V^(r^) is large enough. 

Assuming that the disc Dr is inside the analyticity region of uj{r'^), we decompose 
dDr into 

dDf := dDr n {z : Rez < fi*{lm{z); /3) - 6 Ri{V (T^))} , 

and 

dD^ := dDrn{z : Re2 > iJ,*{lm{z); /3) - eA-^Ri{V{T^))} , 
and write ( |2.2tj| ) as a sum of two integrals If ,„(r^) and /^„(r^) (see figure |I|), 

and 
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2.3.2. An upper bound for /f „(r^). If,^{r^) is not the main contribution to ( p.26| ), 
so that it is sufficient to get an upper bound for this integral Let z G Uq and 
Rez < fi*{lm{z);p). We set 

:= {x E Z'^ : d{x, supp T^) < 1} . 

There exists a constant C5 such that IF^I < CslF^I. iFiom ( p.l7| ) we get 

\uj(r^)\ < exp [ - PWT^W + p\Rez\Ci\r^\ + 3CoS\r^\\ , 

and by the cluster expansion method 

0.(r^)(A(r^)) 



0i(r2)(A) 
We set 

(■= Z- fi . 

Therefore, there exists a constant Cg so that 

|(/'A(r^)| <e-^lir'll(^-^«^"IWp-) if ReC</x*(lm(C);/3)-/i* 
This upper bound imphes 



(2.29) 



j9 ^y"^^ < ^-nf3\\r^\\{l-CeS-rCip-^) 



dD?., 





U*(0;/3) 1 






6 







eA-^Ri{v{r^)) 

Figure 1. The decomposition of the integral into /^^(F^) and /fc„(F^) 



2.3.3. Lower and upper bounds for /^„(F^). In order to apply the stationary phase 
method to evaluate /^„(F^), we first rewrite 0A(r^) in the following form, 

0^(F2)(2) = 0;(f2) e/5Ay(r2)(c+g(r2)(c)) ^ 

where g(F^) is an analytic function of C in a neighborhood of C = and g(F^)(0) = 0. 
Let 

/i*(lm(z);/3) -eA~^Ri{V{T^)) < Rez < ^u* (lm(z); /3) + 6 R2{V {V^)) . 

In this region (see figure [l|) we control the weights of contours with boundary con- 
ditions ip2 and ^1, whose volume is smaller than V^(F^). By the cluster expansion 
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method there exists an analytic function g(r^), which is real on the real axis, so 
that 



exp 



/3n{(pr2 \ip2) + log ^ + log 



e2(intir2) 



ei(r2)(A) 



01(r^)exp (3AViT')C 



d_ 

dz 



=/3Ay(r2)g(r2)(c) 



For large enough /?, t\(3) > T2{/3; 6'), so that we get from Lemma ITT] and Proposi- 

Cnir^i - 



tions ^TT] to pT3 

, d 



dC 



(2.30) 



g(r2)(C)|<2C3i^e--(^^^')(] 
<C7e--(^^^') + C8 



+ 



A Av(r2) Ay(r2 
r^i 



Ciir^l 
Av(r2 



for suitable constants Cy and Cs- Moreover, there exists a constant Cg so that 
(2.31) exp [ - f3\\r'\\{l + Cg^)] < <Pl{r') < exp[-/3||rl(l - Cg6)] . 

Let 

c(n) := npAViT'^) . 
We parametrize (91?^ by z := fi* + re*", — ai < a < 0^2, < ttj < vr. 

/L(r^) - '^^^ 



,c(n)rcosa+c{n)Reg{r2)(C) [ cos(7/;(a))] rfo ^ 



where 

ifj^a) := c(n)r sin a + c{n) Img(r^)((^) — /ca . 
We search for a stationary phase point Cfc,n = ''"fe.ne*"'"'" defined by the equations 

^ fc(n)r cos a + c(n)Reg(r2) (re*")) = and ^V^(a) = 0. 



da V 



These equations are equivalent to the equations ( ' denotes the derivative with 
respect to Q 

c{n) sina(l + Reg(r2)'(C)) + cos aImg(r2)'(C) = 0; 
c(n)r cos a(l + Re g(r2)'(C)) - r sin aImg(r2)'(C) = k. 
Since g(r^) is real on the real axis, ak^n = and r^ ^ is solution of 
(2.32) c(n)r(l+g(r2)'(r)) = A:. 

Lemma 2.3. Let a.i >'it/4:, i = 1,2, A< 1/25 and c(n) > 1. //g(C) is analytic in 
( in the disc {( : |C| < -R}, real on the real axis, and for all ( in that disc 

then there exists ko{A) G N, such that for all integers k, 

k G [koiA),c{n){l-2y/A)R] , 
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there is a unique solution < n < R of ( |2.32| ). Moreover, 

Z ^ < _ / e^(n)rcosa+c(n)Reg{r2) rcos(^(c^))l 



< 



'c{n)rk. 



Proof. Existence and uniqueness of rk,n is a consequence of the monotonicity of 
r I— i> c{n)r{l + g(r^)'(r)). The last part of Lemma 2^ is proven in appendix of 0. 
The computation is relatively long, but standard. 

□ 



Setting c(ra) = n(3^V{V'^) and /2 = 6 /^-^ R2{V {V^)) in Lemma U we get sufficient 
conditions for the existence of a stationary phase point and the following evaluation 
of the integral /^„(r^) by that method. Since rk,n is solution of ( |2.32|) , we have 

kA k , ^ k , kA 



k 



A) 



[I -A) 



and 



cH|g(r2)(r,,„)|=c(n) 
Therefore Lemma ^]3| implies 
(2.33) 



g{r'y{C)dC < Ac{n)rk,n < k: 



A 



A 



'''^c'ic{nr^<pi{Tr<iu^') 



lOVk 



k^ 



^ ^/l + A , , /c! e'^ 



Vk 



2A 



with 

(2.34) c±(A) := (1± A)exp [± 

If A converges to 0, then c± converges to 1. We assume that (see ( p.30|) ) 
(2.35) 



p2| 

A can be chosen as small as we wish, provided that /3 is large enough and -^t^ 
small enough. 



y{r2) 



2.4. Lower bounds on the derivatives of the free energy at finite volume. 



We estimate the derivative of [/Ijlf.^ for large enough k. The main result of this 



subsection is Proposition 

Let < 6* < 1, A < 1/25, and set 



Let e' > and X2 so that 
(2.36) 



e ■=e{l-2^). 

(1+£')X2>X2(00). 
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The whole analysis depends on the parameters 9 and e' . We fix the values of 6', and 
e' by the following conditions, which are needed for the proof of Proposition p73 . 
We choose < < 1/25, 9 and e' so that 

2.37 ed — — - < ^—^ and ^-^ > 1 . 

^ ' 9{l-2^o) d-1 1+e' 1 + e' d-1 

This is possible, since 



> 1 



Indeed, 



n>2 



\^ 1 fl 

^^^(n + l)\\d 

n>l ^ ' 



1 sr^ 1 1 1 

<1- —,+ > -.{-A =e^ 



2d ^ n\ \dJ 2d 



Notice that conditions ( p.3?| ) are still verified with the same values of 9 and e' if 
we replace by < A < Aq. Given 6', the value of 9' is fixed in Proposition p73 . 
/^From now we assume that /3 is so large that all results of subsections |2.1| and |2]2 



are valid. The value of < A < is fixed in the proof of Lemma 2.5 



Given k large enough, there is a natural distinction between contours such that 
9^V{j'^)R2{y{^'^)) < k and those such that 9(]V{r^)R2iV(r^)) > k. For the latter 
we can estimate J^„(F^) by the stationary phase method. We need as a matter of 
fact a finer distinction between contours. We distinguish three classes of contours: 



(1) fc-small contours: 9 (3V (T^) R2{V (T^)) < k; 

(2) fat contours: for r/ > 0, fixed later by (g^, ViV^)^ < ^ l|r^||; 

(3) fc-large and thin contours: 9 /3V (T^) R2{V (T^)) > k, ^(F^)^ > r] \\T^\ 



We make precise the meaning of k large enough. By Lemma |2.2| V ^ VR2{V) is 
increasing in V, and there exists N{x'2) such that 

R,{V)>^ if V>N{x'2). 
yd 

We assume that there is a /c-small contour F^ such that ^(F^) > A^(x2); and that 
the maximal volume of the /c-small contours is so large that remark is valid. We 



also assume (see Lemma 2^) that k > ko{A) and that for a fc-large and thin contour 
(see ( CT) and (p|)) 

C ir^l Cs A 

V(F2) - p^v(r2)i - 2 ' 

so that |g(F2)'| < A, and 

2.38 < — 

pA{l - Ao)vV{r^)-d ~ 10 

are verified. There exists K{A,t],(3) such that if > K{A,r], f]), then k is large 
enough. /^From now on k > K{A,r], (3). 



24 



2.4.1. Contribution to [f1]^i^* from the k-small and fat contours. Let be a /c-small 
contour. Since V i— >• R2(y) is decreasing in V, MA(r^) is analytic in the region 

{z : Rez < /i*(Imz; /?) + ^A-^i^al^^*)} H Uq , 

where V* is the maximal volume of /c-small contours. V* satisfies 

k 



V*~ < 



OPX'2 



Hence 



Since remark 2J1 is valid, we estimate the derivative of Uf^iV"^) by Cauchy's for- 

-s d ^ 2 

mula with a disc centered at //* with radius A~^(Q'x'^ '^^^ P'^k"'^ . There exists a 



constant Cio such that 



(2.39) 



2M(fc) 



< C 



10 



A 



r2;Intr2 90 



(3—{ex'2)— 



k\ /cd-i . 



Let be a fat contour, which is not fc-small. We use in Cauchy's formula a disc 
centered at yU* with radius 

eA-\2{i)v{T^)--^ < eA-^R2{V{T^)). 

We get (see ( ITTUD ) 



A(rY](;i|<fc! 



< k\ 



( Ay(r2)f 
X2{i)e 



-n[ri(/3;0')-C5<5]|r2| 



p2|^g-n[ri(/3;e')-C55]|r2l _ 



We sum over n and over using the inequality 

Y^mJ>e~'''^<-J{p + l) (p>2,g>2). 



m>l 



There exist Cn and Ci2(^') > so that 



r^iintr^'so 
V(r2)^<,,||r2|| 

F'^ not fc-small 



We choose r] so small that (see (^ 
(2.40) 



A(C2r?)^ 
.(Ci2/9)^X2(1)^, 

A(C2^)^ 
.(Ci2/5)^X2(1)^. 



A;!r 



c/- 1 



A 



< 



< 



A;! k— . 



A 



A(C2^)~ 

(Cl2/3)^X2(1)^ /?^(^X2(00))^ /3^(^X'2)^ 
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2.4.2. Contribution to [f1]^i^* from the k-large and thin contours. For fc-large and 
thin contours we get lower and upper bounds for [(j)\{T'^)'^]^l^}. There are two cases. 

A. Assume that Ri{V{T^)) > R2iV{T^)), or that V{T^) is so large that 

ef]v{r^)Ri{v{r^)) >k. 

For each n > 1 let c(n) = n/SAVlT"^). Under these conditions we can apply Lemma 
^ with a disc „ so that (9D,, „ = dD^^^. Indeed, if Ri{V{T^)) > R2{V{T^)), 
then we apply Lemma with R = 6A^^R2(y(r'^)), and in the other case we 
set R = 9A~^Ri{V{T^)). In both cases rk,n < R, which implies 9A-,,„ = dD^^^. 
Therefore we get for /^„(r^) the lower and upper bounds ( p.33|) . 

Lemma 2.4. There exists a function D{k), Yrnvk^^Dik) = 0, such that for j3 
sufficiently large and A sufficiently small the following holds. If k > K{A, t], (3) and 
Ri{V{T^)) > R2{V{T^)) orepV{T^)Ri{V{T^)) > k, then 

(1 - D{k)) [0A(r2)]J2 < -MT')fJ < (1 + D{k)) [0A(r2)]Ji . 

Proof. We have 

n>2 [YAyi- )\fj,* 

From ( |2.33| ) there exists a constant C13, 

^2\n-\{k) 



The isoperimetric inequality (|2.21| ), R2in) < X2{n')n~^ and the definition of /c-large 
volume contour imply 

^lir'll > PX2{V{T^))V{T^)'^ > epR2{V{T^))V{T^) > k. 

Let b := Cg6 (see ( p.31| )); we may assume jQ—b> | by taking P large enough. Then 

~ n>2 ~ n>2 



— k / J 

- n>2 

< (^e-l^)'5^e->^ 



n>l 



We choose A so small that c+{A)c-.{A) ^e «> < 1. 



B. The second case is when 

e(3V(r^)Ri{v(r^)) <k< e(3V(r^)R2{v{r^)) . 



□ 
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Since the contours are also thin, 

(3\\T'\\<r'o-\,{i)-'pexi{i)v{T'Y-^ 

< r'e-\,{ir'(3ev{T^)R,{v{T^)) 

< r]~^e~\i{l)~^k = Xk. 

We choose R = l3A~^R2{V{r'^)) in Lemma The integration in ( |2.26D is de- 
composed into two parts (see figure |I]). We show that the contribution from the 
integration over dD^^ is neghgible for large enough jS. Since k > i^(y4, 77,/?) and 

the contours verify V{T'^)^ > ?7||r^||, we have 



A{1- A)7]V{T^)^ A(l - Ao)r/\/(r2)d 
By definition of K{A,ri,p) (see (|23|)) 

n(3\\T'\\p"'Cin,n<^. 



From (|2.29|) with r = rfc .„ we obtain that the contribution to |[nA(r'')]^^i| is at most 
(1 + (/?Ay(r^))%xp (A) 1^ Y: n^e-^ll^^ll(-^«^) . 

n>l 

As in the proof of Lemma |2.4| , we choose jS large enough so that we can assume that 
^-C<,6> l Then 

n>2 

<g-/3||r2||(l-C65)j^l + ^g-i 



n>l n>2 



10"'^ 



n>l 

= e-/3||r^ll(i-C65)(i + ^(^))_ 

Since /3||r^|| < A A;, by choosing A small enough and /5 large enough, so that 6 is 
small enough, we have 

(1 - D(A;))c^e'=e-^ll^'ll^«'^ > (1 - D{k))cte''e-''^^''^ > ef 

and 

(1 + D{k)){l + A)'=e^e^ll^'ll^'"^ < (1 + D{k)){l + A)*^eise^'=^«^ < et . 

If these inequalities are verified, then the contribution to —[uxiT'^)]^!^} coming from 
the integrations over dD^^ ^ is negligible with respect to that coming from the inte- 
grations over dD"}^ ^ . Taking into account ( |2.33| ) we get Lemma 



Lemma 2.5. There exists < A' < Aq so that for all (3 sufficiently large, the 
following holds. If k > K{A',ri,P) and is a k-large and thin contour, then 

-K(r^)]?*^ > ^{^-D{k)){f3AV{T')fctct>l{T'). 



27 

Proposition 2.4. There exists (3' so that for all (3 > (3' , the following holds. There 
exists an increasing diverging sequence {A;„} such that for each kn there exists A(L„) 
such that for all A D A(L„) 

Ci4 > is a constant independent of (3, kn and A. 

Proof. We compare the contribution of tlie small and fat contours witli tliat of the 
large and thin contours for k > K[A' ,r], f5). The contribution of the small contours 

to |[/a]^*| is at most 



Cio (3-— {ex2y~ k\ k— < Cio A'' (3-— 



ed 



Ox'2 



"d-l d 



The contribution of the fat contours is much smaller by our choice of t] (see (|2.40 )). 
The contribution to — [/a]^;^* of each large and thin contour is nonnegative. By 
assumption ( |2.36D and the definition of the isoperimetric constant X2, there exists a 
sequence F^, n > 1, such that 



lim lir^ll ^ oo and ^(r^)— > 

n—^co 

Since x ^-le"^ has its maximum at x = k-j^i^ we set 

d-l 



ir? 



kn 



d 



For any n, is a thin and fc„-large volume contour, since by ( |2.37| ) 



(3 {I- 2^')V{T^)R2{V{T^)) >(3{l- 2y/A')V{T 



2\ I 



If A D F2 then 



^ l-D{k) ^_±^ 

> )-J-A''"B ~ 

-20 ^ 



dc_ 



{d-l){l + e')x'2 



and (see ( pIMD ) 



kt'M^l) > kt' exp [ - [K j—^ + 1) (1 + 0,5)] 

By the choice (p.37|) of the parameters 9 and e', if 5 is small enough, i.e. /3 large 
enough, then 

1 

d c_ 



1 

ed 



d-l 
d 



< 



-C95 



e{l-2^/A') d-ll + e' 
Hence the contributions of the small and fat contours are negligible for large kn (see 
( p.39| ) and ( p.40|) ). Let A(L„) be a box which contains at least |A(L„)|/4 translates 
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of r^. For any A D A(L„), if kn and /5 are large enough, then there exists a constant 
Ci4 > 0, independent of /3, kn and A D A(L„), such that 

□ 

2.5. Lower bounds of the derivatives of the free energy at infinite volume. 

We show that we can interchange the thermodynamic limit and the operation of 
taking the derivatives, and that the Taylor series, which exists, has a radius of 
convergence equal to 0. These statements are a consequence of Lemmas |2.6| and | 



Lemma 2.6. If j3 is sufficiently large, then for any k there exists Mk = Mk{(3) < 
oo, such that for all t G {fi* — e, fi*] and for all finite A, 

\[fl]['^\<M,. 

Proof, ujir"^) is analytic and ri(/3, ^')-stable on a disc of radius ^A~^i?2(^(r^))- 
From Cauchy formula 

|[nA(r2)]f I < A;!C^5|r2|3^e-^'^|r'l, 
for some constants C15 and k > 0. Therefore 

r^cA r^cA 



□ 



Lemma 2.7. 



hm [/l(,)]^i = lim[/]('^) 



Proof. We compute the first derivative at the origin. Let rj > 0. 

/(/i*)-/(/i*-r/) 



V 

L^oo rj 

j.^ [/a(L)] V + MfAiL)fJ-xUv) 

By Lemma l[/A(L)]|f*-a;i(r;)l — Therefore {A{ri)}n is a Cauchy sequence. 

Hence the following limits exist, 

1/1';.' = lim ^'''''- f^'-"> = li„;[/l« = . 

Same proof for the derivatives of any order. 

□ 
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